* + x 


Poderíamos também considerar um elemento que embora transmitindo uniformemente todas 
as frequências no intervalo [0, F] introduza um atrazo «T». 
O operador linear de transmissão relativo a este quadripolo será: 


R ()= E qd W<o<W (39) 
R ()=0 Io |>DW 
Obtém-se para este caso 
Fa (t) = PESE, Ses sen W (t—r) (40) 
= (t—s) 
A série uq reduz-se também a 
Uç — 
u4=0 q5o 


desde que se arbitre To = 7 


8.1.4 — Elemento inverso 


Elemento inverso dum dado sistema será, por definição, um quadripolo que associado em 

cadeia com o referido sistema conduza a um quadripolo equivalente que é o elemento unidade. 
Dado um sistema caracterizado por [un, 71] definamos o elemento inverso por [va, 72). 
Determinemos pelas formúlas (32) a série [u.p, To ] equivalente ao conjunto. 


Teremos imediatamente: 
co 


gm & VYo=a ta | (41 a) 
n=— 00 


w=u+m7 (41 b) 


Como uep se conhece, pois trata-se da série correspondente ao elemento unidade, (41 a) traduz 
um sistema de equações que nos permite determinar vg == vp-n. Embora se trate dum sistema com 
um número infinito de incógnitas práticamente reduz-se a número finito. 

As equações (41 a) e (41 b) reduzem-se a uma forma simples desde que se escolha T1+73=0 ==, 
pois que, nessas condições, é como se viu 


Uco == 1 
Uep == O pão 
Virá então: 
ps 
à Wentn=1 p=0 (42) 
n=— co 
rs 
& VYp-nh=0 po 
n=— 00 


8.1.5 — Redução de um quadripolo qualquer a um quadripolo de característica A (5) 
limitada a [0,F]. 


Embora se trabalhe com funções de espectro limitado acontece que por vezes os quadripolos 
não têm característica de amplitude limitadas a [0,F]. Nessas condições já a resposta ru (t) não é fun- 
ção de espectro limitado, invalidando os resultados obtidos. 
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Podemos facilmente reduzir esse quadripolo a um quadripolo com A (m») limitada a [0,F]. 
Bastará associa-lo em cadeia com o elemento unidade. 

Dessa associação resultará, como é evidente, um quadripolo que no intervalo [0,F] tem um 
operador linear de transmissão que coincide com o do quadripolo primitivo e que se reduz a zero 
para as frequências superiores a F. 

O facto resulta do elemento unidade funcionar como um filtro passa-baixo ou ainda da consi- 
deração de que o operador linear de transmissão do quadripolo equivalente à associação em cadeia 
é o producto dos operadores de cada um dos elementos. 

O problema que se põe será pois de determinar a série us correspondente à resposta ru (t) 
filtrada pelo elemento unidade. 

Um processo será, evidentemente, determinar em vez da resposta 


(= [ DR6) cido (43) 


a resposta filtrada ru (t) que corresponde a integrar-se entre [—- W, W] 


a W — e It 
EE (t) = | im R (0) é do (44) 
U -W 


Outro processo será determinar a resposta que o filtro passa-baixo [0, F] dá à função ru (t). 

Basta aplicar a expressão que nos dá sz (t) em função da resposta r, (t) que o elemento uni- 
dade dá ao impulso unitário. 

Usaremos a expressão considerada no capitulo 3. 


sa (t) =[" se (ES) ra (1) dr (45) 


Simplesmente aqui, a integração deve estender-se de— co a + o pois que sendo o elemento 
unidade fisicamente irrealizável, não satisfaz a 


Tu (T)=0 para 1 < O 
existindo ru (7º) desde —- co à co, 


Pelo resultado obtido em (35) r, (t) é dado por 


1 1 
r, (t)=——senWt=-—sen 27 Ft 
Tt nt 


Nessas condições a resposta filtrada será dada por 


"o 1 
ro (t) = | Tu (t—7) —— sen 27FErdr (46) 


“ — co 


Qualquer dos processos conduz à função r! (t) que é uma função de espectro limitado [0, F] 
e à qual corresponde uma série u, que permite aplicar ao quadripolo considerado a técnica das 
séries no tempo. 

Indicaremos, quando tratarmos das aplicações, um modo prático e simples de obter u, sem 
sem ter que se recorrer a (41) ou (46). 
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8.1.6 — Funções de espectro não limitado 


Por vezes a função si (t) não é de espectro limitado mas destina-se a ser transmitida por um 
canal cujo operador linear de transmissão é limitado a [0,F]. 

Nessas condições deve começar-se por reduzir si (t) a uma função sr (t) para o que basta til- 
trá-la através do elemento unidade. 

Aplicando (46) obtém-se imediatamente a função si (t) filtrada, sf (t), 


sE (t) = oi S1 (t— 2) =—sen 27FET dT (47) 


— (O 


Como o canal está limitado ao intervalo [0, F] a resposta que se obtem é a mesma quer se 
considere s; (t) ou st (t). 
À função de si (t) pode, agora, aplicar-se a técnica das séries no tempo. 


10 — RESUMO DE RESULTADOS. EXPRESSÃO SIMBÓLICA DAS SÉRIES NO TEMPO 


Os conceitos introduzidos e resultados obtidos constituem os elementos fundamentais da teoria 
das séries no tempo relativa à sua aplicação à determinação da resposta de sistemas lineares de 
transmissão. 

A reunião dos referidos elementos constitui, pois, a referida teoria. Antes de prosseguir 
entrando no domínio das aplicações faremos um resumo dos resultados obtidos. 


10.1 — Condições de aplicação 


A teoria estabelecida apenas se aplica a funções, sr (t) de espectro limitado ao intervalo [0, FF] 
e a quadripolos cujo operador linear de transmissão seja limitado ao mesmo intervalo. 


10.2 — Séries no Tempo. Definição 


São as séries cronológicas cujos termos são os valores da função sr (t) tomados nos instantes 


«tn». Isto é 
Série no tempo == sr (tn) 


tn == nTr + T Tr=— 
7 — denomina-se tempo de referência. 
10.3 — Séries normalizadas 


São as séries no tempo em que os termos são os elementos sr (tn) divididos por 2 F 


= SF (tn) 
ns a 
10.4 — Propriedades fundamentais das funções sr (t) 


Uma função de espectro limitado, sr (t), e o seu espectro, E [sr (t)], exprimem-se completa- 
mente em função das correspondente séries, sr (tn) OU Sn 


sao FE E MR pres 2 m —— sm 
n=— 0 2%EF (t—t,) n=— co 7 (t—tn) 
is 1 pad — jwtn Eis — jutn 
Propriedade II E [ss(U]=—  Z sr(tn): aÃ o ta é 
2W n=—co 2% n=—o0 
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10.5 — Princípio fundamental 


Um quadripolo em que o operador linear de tronsmissão esteja limitado a [0, F], apresenta uma 
resposta ru (t) ao impulso unitário sy (t) que é uma função de espectro limitado [0, F). O quadripolo 
fica completamente definido pela série correspondente a ru (t), designando-se por un a série 


normalizada 
o Tu (tn) 


2F 


+ + + 


A resposta sz (t) a uma função de entrada sy (t) obtém-se então a partir das séries correspondentes 


as, (t)eru (t). 
A série normalizada sz p correspondente a sa (t) obtém-se a partir de (31 a) e (31b) 


o 


Sap — 2 Up — | Sin (31 a) 
n=— 


Ta == To + KS! (31 b) 


Designando simbólicamente 
a série sn por 5 


e convencionando definir (31 a) e (31 b) como o produto de U por Si podemos escrever simbolicamente 


So == U.S: (48) 
* * 


Dados dois quadripolos de séries un e u2q referentes às respostas rui (t) e ruz (t) podemos 
escrever, dentro da mesma simbologia, para a série equivalente U 


U = Ui. Us (49) 


A série equivalente é o producto das séries de cada elemento, dado que up se obtém de un 
e u:y por expressões (32) formalmente iguais a (31). 


+ 4 + 


Estas relações traduzem o princípio fundamental que preside à aplicação das séries no tempo 
ao problema proposto, determinação da resposta de sistemas lineares de transmissão. 


11 — TEORIA ELEMENTAR DAS SÉRIES NO TEMPO 


Podemos, agora, com os elementos de que dispomos reduzir a teoria feita a um conjunto de 
propriedades que convém reunir de modo lógico e sucinto. 
Teremos assim: 


(11.1) Notações 


Si — expressão simbólica de sm, série normalizada de sy (t) 
S: expressão simbólica de szp, série normalizada de ss (t) 
U expressão simbólica de u,, série normalizada de ru (t) 
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1 expressão simbólica da série normalizada da resposta ru (t) do elemento unidade. 

U-1 expressão simbólica da série normalizada da resposta ru (t) do elemento inverso. A série 
va correspondente ao elemento inverso determina-se a partir de (41) 

Si expressão simbólica da série normalizada correspondente à função si (t) resultado de 
filtrar pelo elemento unidade a função si (t) de espectro não limitado 

U* expressão simbólica da série normalizada correspondente à resposta ri (t) resultante de 
filtrar pelo elemento unidade a resposta ru (t) dum quadripolo de operador não limitado 
ao impulso unitário iu (t). 


11.2 — Propriedades fundamentais 


Definido o produto de séries pelas expressões (31), teremos fundamentalmente as seguintes 
propriedades das séries no tempo cujo conjunto constitui a teoria elementar das mesmas: 


I — Lei de transmissão de sinais 

Sa = U.S: 
IH — Lei de associação de quadripolos em cadeia 

U=U,.U; 
HI — Elemento unidade 1 

U.1i==U 
IV — Elemento inverso U- 

EA Sir E | 


V — Transmissão de funções de espectro não limitado através de canais limitado a [0, F]. 
Usa-se a série Si 

VI — Transmissão de funções de espectro limitado através dum canal não limitado a [0, E]. 
Usa-se a série U* 


Damos, aqui, por estabelecida a teoria das séries no tempo relativa à aplicação que encontram 
na determinação da resposta de sistemas lineares de transmissão. 

Podem-se ainda estabelecer outras propriedades e resultados que se deduzem lógicamente dos 
princípios estabelecidos pelo que incluiremos estes últimos no domínio das aplicações. 


(Continuação) 


TECNICA 


C. D. U. 624.044.2 


O METODO DE RELAXAÇÃO 


SUA APLICAÇÃO AO CÁLCULO DE ESTRUTURAS RECTICULARES 
CONTÍNUAS PLANAS CONTENDO PEÇAS CURVAS 


(Conclusão) 


PELO ENG.º ClVIL ANTÓNIO PACHECO CORREIA DE CINTRA 


II APLICAÇÃO 


Forma, dimensões e solicitação 


Elementos o Secção | E | ! 
Pilares | Arco (m?) | (mt) (Tn/m?) (m) 
| 
1,2 — 0,30 >< 0,30 > 10 | Mm o 4,00 
3,4 - 0,30 X 0,30 6,7 >< Lo | 21x 105 | 4,00 
- 2,4 | 0,30 X 0,30 É,7:>< 16 21 >< 10" 15,00 


Condições de apoio 


Encastramento perfeito de ambos os montantes. O montante B sofre um assentamento desi- 
gual de 0,005 m, 
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as 
nm + 

a 
A = 
A bor ENT 


O MATERIAL DE va Lmestent ACEC 


CIENTIFICAMENTE ESTUDADOS 
E VERIFICADO, OFERECE AS MAIS . 
COMPLETAS GARANTIAS = ALTA 


Q U A .: D A D É 


E. “LAMPADAS 
Wo e BALASTROS 


ES + 


A 
Sa 

cs R 

z . te “> 
Poe 
“a 


UMA DAS FASES DE FABRICAÇÃO 
DE LÂMPADAS FLUORESCENTES 


ACECÇ 


ACEC - ATELIERS DE CONSTRUCTIONS E a - 
ELECTRIQUES DE CHARLERO! Eae qa al | 
REPRESENTANTES GERAIS E CHARLEROI: 


INEL - INDÚSTRIAS ELÉGTRICAS ASSOCIADAS 


SEDE EM LISBOA: RUA RODRIGO DA FONSECA, 76, 4.º 


DELEGAÇÃO NO PORTO: RUA DE SANTA CATARINA, 470 
TELEF. 28841, 28842 e 28843 — TELEG.: INELAPO-PORTO 


TÉCNICA — XXI 


a COMPANHIA UNIÃO FABRIL 


põe à disposição dos seus clientes 


no Estaleiro Naval 


Rae pa e4] Sá 


Geigre 


— Uma oficina de Ensaios e Reparação Diesel 


— Uma máquina electrodinâmica de equilibragem 
para rotores de 30 a 3000 Kgs. — 2 M. de 


diâmetro máximo — 3,420 m. entre apoios. 


ROCHA DO CONDE DE ÓBIDOS 
LISBOA Telef. 6621 48 


TÉCNICA — XXI 


o=+>xg Ega —g'LeS + 1't6z + SrEz: —(S) 

o = bx g'Egz — g'lzS + S'têz + 1'E6c: — (b) 

o=+t>g'lzS + SrxS'L€ +SEoL—ziSzc—higã: — (€) 
o=+> g'LzS + SrxS'LE + SEolL — b'iga + ciScc—:— (2) 
o=oS>Sr —-SLEACLE:— (1) 


ogdpIGdO Op 8910/92] O BIQUGN|UI OP Sojua]0909 SOp OBÍBIHIJOA 
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oCL |oSL—| g'lzS — | g'lzS — | SEolL — | SEolL — || g'laS g'LeS 96991 — 9'696'1 — 1 == 14 = és 9 
o o g'Ega g Egg g LzS g'lzS gfgz - | gfgz— | gLlzs g'LiS I='N=n L 
o o g'lzS — | gLazS SfoL — | S'€oL g'gs 9 88 — g'ietSc — | gatSc [I=uU=tH— 9 
o o g Egz o g'LaS o1gê — | c'li1 1'E6z Ctêz I=9% S 
o o o ide o vas | c'lii o IgE — Stêz 1't6z 1 = 81) + 

oLE |SLE—| g'laS — SEL —| o. | rE6z Gh£z ciSec— | bige p= E 

oLE |SLE — o 9' les — o SfolL — | S+êz 1'E6z b'iga ciiSee — 1 = by z 

oS — lo'S o a o o o o S'LE S'LE [1 = 84 I 
E ó end do EN IN g “X 'N EN 

odiL eN 
oI 6 9 S y (3 c I 
ogôeiado op S910/9E4 

tuge | =btg S'LE S'hêz o'S eli | ooSr | orxL'9| OLX IzZ b'z 
S £oL o'Lot'1 - g'LeS - g'Egz 00 + Ot L'Q| OLX IE he 

S £oL o'Lob:i - g'LeS - g'fgz | oo'y ;-01 X L'9 | OLX IZ | E 

J us [JÁ | [38] mi qu) | (um) — Guug) | 
-— -—— | “Wa 
É ba a DE etica dO CU ; 


— ——— ————— ————— —————e >> >>> >> 


BIQUGN|UI Op S9juo/9209 


Reacções nos extremos dos elementos considerados perfeitamente encastrados 


gl” 1x 16 


Mi=M=-'—= ——— = 1,333 Im. 
: 12 12 am 
Montante 1,2 
| 15€4 
Hi me Ho — q .= ——————. == 2 Tn. 
2 2 


Fig. 17 


Acções desiquilibradas dos nós 


Fig. 18 
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Relaxação 


6 000165 |. 


Valor 
Tipo da 
“ relaxação 
I o 00500 
7 0,00758 
2 0,00312 
> | "Q0UGO7 


2 — O 00031 


3 — 0 00005 


0,00091 


8 0,00012 


7 0,00026 


3 O 00007 


2 0,00007 


3 0,00001 


5 0,0001393 


gm 0,00004 


| 0,00003 

Di | 0,00003 
7 | 0,00007 

6 | — 0,00002 

3 | O 00004 


O 000TO 


n93% 


— 7024 


o,co8 


— O 720 


— 0,712 
0,698 


— 0014 


— 0,028 
0.267 
0,239 


— O 296 


0,003 


0,137 


O,I40 


0,020 
0,160 
— 0,158 


0,009 
0,003 
0,005 


0,029 


0 034 
9.030 


0,064 


0,012 


WTC TT —ee—e Ni 


Ny o) 
0,0CO 2,000 
0,187 0,000 
0.187 | 2000 - 
| 3999 | —2000 
| 4186 | 0,000 
Din. Zi A hm 7 DO 
O 007 — 0,097 
| 0,878 o,914 
0885 | 0,817 
— 0.900 — 0,360 
—ooI15 |—0457 | 
'— 0,087 | —oogr 
“— o, 102 0,366 
O,T13 — o,012 
OO! | 0,354 
| 0213 | —0347 
0,224 0 007 
— 0236 | 0,063 
“— 0012 | 0,070 
0,137 | —0069 | 
0,125 0,001 
— 0.158 o,016 
— 0,033 0017 
0.020 o 020 
= O. | 0097 . 
— 0,025 | 0.002 
— 0,038 | 0,039 
0,023 — 0,038 
ema oor5 | 0001 | 
0.038 0015 
— 0,023 — 0016 
0.009 — 0,015 
0,092 | ooo | 
0,008 0,009 
0,040 | oolo 
— 0 067 O 007 
— 0.027 0,017 
0 037 | -— GIO: 
0,010 — 0,001 
0.051 | o 00! 
0,061 | 0000 | 
— 0 090 o 009 
— 0 029 0,009 


RS E-S 


— 2,000 
— 2000 

0,097 

— 1,903 

0,732 
TIE | 

1,170 


— 0,001 


Verificação da relaxação 


| Valor I | 2 3 4 
Tipo | da | | . 
| relaxação No Ni X» X4 
888. | o oco 2,000 oO 000 
o 00500 0,187 0,187 0,000 0,000 
2 | 0002! | —S5st | oBig | oBs2 | o68 
| 3 | * 000010 | o 028 — 0.227 | 0,022 0,029 
“4 | ocoos || 0.308 0,245 — 0,400 | O,124 
E | — 000294 | — 0,690 — 0862 | — 0,345 I,120 
E 20 0,00163 | 4,128 — 4,128 | —o,096 o 096 
2 | ooo | qu | 4173 | — 2087 [208] 
8 |. oGsia | — 0,236 | —o236 | o 0,063 0,063 . 
E — gom — 0 029 | O 009 O 027 


Valores totais das rotações e translacções 


rz == 0,00140 rad. 
r4 == 0,00185 rad. 
Vp == 0,00500 m 
Us == 0,00895 m 
U, = 0,00496 m 


Verificação das reacções de apoio 


Dos valores obtidos na fig. 19 resulta 


| 0,096 — 0,096 = O 
Equilíbrio de forças | 
4 — 3,624 — 0,335 —- 0,041 


Equilíbrio de momentos | 4><2 — 0,096 x 15 — 5,075 — 1,321 == 0,064 


Deverá ainda proceder-se à verificação dos equilíbrios de forças e momentos nos nós. 
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Determinação dos momentos e esforços nos extremos das peças 


E 


Ei | [ex] am E) ss | =D) lry las |yy lap (Vo—VA)) xrlas Lyrlas | ÍA E TA Fp (1.3) (2.93) (4.6) (5.6) (7. I1) (7.12) (8.11) (8.12) (9.11) (9 . 12) (10.11) | (10. 12) 
ou al 
A,B (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (2x) (22) (23) (24) 
1,2 | 527,6 | 2638 | 000896) — ne o 5276 | — | 14070 | 7035 o 0,00140 | 4,727] 2363] — a o 0,739 Ra a o 1,970 o 0,985 
34 | 527,6 | 263,8 0,00497 | — -— o 527,6 — T.407,0 | 703,5 o 0,00185 2,622| 1,311 — -— o 0,976 =— — o 2,603 o 1,301 


2, 234,5 | 117,2 | —- 0,00399| 37,5 5,0 0,0050 | 234,5 | 37,5 844,2 | 281,4 | o00or4o | 0,00185 | — 0,935| — 0,467 0,181 0,025 0,328 | 0,434 | 0,052 0,069 1,182 1,562 1,394 0,520 


DE O O O O O O 


À À 


Momentos Esforços nos apoios 
Elem. Mm Pilares | Arco — Pilares ro 
MA Mp Ps Ma Mp o Ka Xa Ya 
AB |—1ig+214-24 |—1g4-22423 | —I3— I5S+21— 24 13— IS +22 — 23 — 14+17+18 | — 14 + 17 — 18 — 16 + 19 +- 20 
1,2 — 3,742 — 2,757 - ss — 1,624 ei is 
314 — 1,321 — 0,019 — — — 0,335 — — 


(qm) (b) (c) 


Fig. 19 


III APLICAÇÃO 


Forma, dimensões e solicitação 


Fig. 20 


Dada a simetria da solicitação e da estrutura, considerar-se-á apenas metade desta no estudo 


que se segue. 


Elementos Secção I E ! 

pil are [e (m*) (mt) CEn/mº) (m) 

1,2 ss — 0,30 x 0,30 | 67x 10"*| 21x I0º 4,00 
3.4 e — 0,30 0,30 6,7x 107* 21x 0º 4,00 
56 — em 0,30> 0,30 67x 1071 21 X 10º 4.00 
== 2,4 mem 0 30x0 30 6,7x 10-* 21x 10º 15,00 
sea 46 — 0,30 0,30 67x 10'*| 21x Io 15 00 
-— — 2.4 0,0002 | -— 21x 10º 15.00 
q — 4,6 0,0002 — 21>< 10º 15,00 


Condições de apoio 


Não há qualquer assentamento dos apoios. O seu encastramento é elástico. Recorda-se, que 
por encastramento elástico se entende o encastramento, em que as rotaçães e os momentos reacção 
que se desenvolvem são proporcionais. Assim, se conhecermos o valor do momento que se desen- 
volve M,, para uma rotação unitária, teremos possibilidade de conheces o momento desenvolvido, 


para qualquer rotação diferente da unidade (1). Seja o nosso caso: 


M, = 250 Tm 


Valor de R« para os tirantes 


Será como na aplicação I: 
Ri=280 Tn 


(1) — Sobre este assunto consulte-se «Cálculo de Estruturas Atendendo ao Grau de Encastramento das Funda- 
ções», pelo Eng.º Brazão Farinha. «Boletim da Ordem dos Engenheiros» Vol. V, n.º 8. 
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Coeficientes de influência 


un sia 
Elem E I ! = io | rx | lry| rr 
[xx | yy N peste a 
A,B (Tn/m?) (mt) (m) xr | | |yr| m f 
12 | 21x 10º |67x10-*| 400 | 263,8 o 527.6 — | — 1.407,0 
24 | 21><10º |6,7>x107*| 15,00 117,2 50 | 234.5 Edy | 844,2 281,4 
Factores de operação 
I 2 3 4 5 | 6 7 | 8 
N Tipo | | 
N No | XK | % No | Xj Va | 
I n=1|—1.657,2| —70o35| 5276 |— 527,6 0.0 0,0 00 | 00 
2 =2) —7035|— 2.251,2| 2931 | — 527 6| 281,4 |, 2345 |—3%5 375 
3 Us= 1 527,6 293,1 |— 661,0 | 2638| 2345 397,1 | 0,0 00 


Verificação dos coeficientes de influência e factores de operação 


(1) —: — 1.651,00 -- 703,5 + 527,60 4 — 2500 — O 
(2)-—-: — 703,5 — 2.251,24 281,4 + 4x 227,60 + 375xX I5=0 
(3) —: — 527,6 + 293,1 + 234,5 — 2638 x 4 = 


Reacções nos extremos dos elementos considerados perfeitamente encastrados 


Me = Ma mens O 
|2 3x 225 
E E es = Ely: o: E = 28,125 Tnm. 
Arcos 24 e 4,6 | E 8 x 
| 3> 15 
Ve = Vi-T—-=————— == 22,500 Tn. 
| 2 


9009 
28,125 


Fig. 21 
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Acções desequilibradas nos nós 


22,500 


28,125 28125 TN 
2 19000 9000) 4 


1 3 


Fig. 22 | 
Relaxação 
Valor I 2 3 
Tipo da 
rélaxação N N; Ko, 

0,000 0,000 — 28,125 

3 — 004254 | — 22444 — 12,468 — 28,119 
oo — 22444 | — 12468 | — 0,006 . 

I — 0,01354 22,435 9,525 — 7144 

a far”) — 0,008 — 2,943 -— Qui 

2 — 0,00130 | ogI4 2,927 — 0,381 

daiadd 0,906 — 0,016 — 7.531 

3 — 0,01139 — 6,009 — 3.338 7,529 

— 5,103 — 3.354 — 0,002 

I — 0,00307 5,087 2,160 — 1,620 
PRO À ce DO | — LIM | — gm 
2 — 0,00053 O 373 1193 | — oI55 

0,357 -— GOA ans, 

3 — 0,00268 — 414 — 0,785 1,771 
Dolo | — os | — 07386 | "0006. 

Ee) — 0,00063 LO44 | 0,443 — 0,332 

| — 0,013 — 0,343 — 0338 

2 — 0,00015 O,105 o 338 — 0044 

0,092 -— 0,005 -= 0,382 

3 — 0,00058 | — 0,306 — 0,170 0,383 

| — 0,214 — 0,175 0,001 

I — 000012 0,199 o 084 — 0,063 

— O,0I5 — 0,091 — 0,669 

2 — 0,00004 0,028 0,090 — OoIz 

0,013 — 0,001 — 0,074 

3 — 0,00011 — 0,058 — 0,032 0,073 

— 0,045 mu. QUI — 0,001 

I — 0,00002 o 033 0,014 — OoIo 
CC [lo | — co2 | — oog | — oo. 

2 — 0,0920008 0,006 0018 — 002 

— 0006 | — 0,001 — 0,013 

3 — 0,0000077 — 0003 — 0,002 o 005 


Verificação da relaxação 


a RasoE K à 3 
ii | relisação Ni No X> 

| 0,000 0,000 | — 28,125 

I — 0,01738 28,799 | 12226 | — 9,169 

2 — 0 002028 1.426 4 566 — 0594 

3 | — 0,057307 | — 30234 | — 16 795 * 37880 

— 0,009 | — 0,003 | — 0,008 


Valores totais das rotações e translacções 


ri==— 0,01738 rad. 
ra == — 0,002028 rad. 
Us == — 0,057307 m 


Exforço no tirante 


AIE. >< 21x 10" 
pm DORES o DONAS RI gta 


Verificação das reacções de apoio 


Dos valores vbtidos na fig. 23 resulta 


Equilíbrio de forças Ei sit: 
É 68 | 3><30 — 22'576 — 44,858 22,576 =0 
Equilíbrio de moment. — 90X 15 + 4,354 — 4,354 + 22,576 X< 30 + 44,858 x 15 = 0,027 


Deverá ainda fazer-se como já temos referido, a verificação do equilíbrio de forças e momen- 
tos nos nós. 
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HEENAN & FROUDE LTD. ——— 


WORCESTER 
ENGLAND 


VELOCIDADE VARIAVEL COM 


Motor «shunt» de corrente continua 
Unidade electrónica de rectificação 


Regulador de velocidade 


Representada por: 


THE ENGINEERING COMPANY OF PORTUGAL, LDA. 
Rua dos Remolares, 12-1.º — LISBOA 


TÉCNICA — XXI 


Tele l 


Caracteristicas: 


Leveza; grande resistência ao choque, esmagamento e abrasão; 
boa flexibilidade; grande estabilidade dimensional, perdas de carga 
minimas;, imunidade à incrustação, insensibilidade às baixas tempe- 
raturas;, resistência aos ácidos minerais, alcalis e soluções salinas, 
inatacabilidade por fungos ou roedores. Não tóxico. 


Restrições : 
A temperatura máxima de serviço continuo é de 50º (, 
Aplicações: 


Construção Civil, Indústria Química, Condutas desmontáveis ou 
tixas para cerveja, agricultura, vinicultura, sulfato para as vinhas, 
estábulos, nitreiras, centrais leiteiras, sumos de frutos, águas ácidas, água 
salgada ou calcárea, condutas de água e ar nas minas, pedreiras, etc. 


SOCIEDADE FABRIL DE MATÉRIAS PLÁSTICAS 


PORTO LISBOA 

RUA DO HEROÍSMO, 2091 RUA DA EMENDA, ' 
gramas UNISOTRA , Tele | gramas UNISOTRA 
fone (prov.) 5 2102 fone 2 0448 - 56 T4BB - 56 7489 


TECNICA — XXIV, 


Determinação dos momentos e esforços nos extremos das peças 


Elem. Txlas | XE | id (Us — UA) | ES AB loé| > XT | AB | JT AB m|FE|AB elrr | AB TA Tg (1,4) (2.93) (4.6) |(6::6) | (7.11) (7. 12) (8.11) |(8.12)| (9.11) (9.12) | (ro. II) | (IO. I2) 
ou AÍ = 
> 
A,B (1) (2) (3) (4) (5) |46)) (7) (8) (9) (10) (12) (12) (13) (14) (15) | (16) (17) (18) (19) (20) | (ax) | em 4 (23) (24) | 
E ed! DAT ED D |———— | — | ERR BU pa af 
1,2 | 527,6 | 2731 | —0,057307 | — = || 8276 | 1.407,0 | 703,5 | —0,01735 | — 0,002028 | — 30,235 — 15,117] — | — —9170 =1070) —— | — |=24454]=2/853| 12,227 | — 1,427 
844,2 281,4 | — 0,02028 o | 13.438 6,716, o o —0,478 | O —0,076] O Pa [o |cogn | O 


sm, 
24 | 2345 | 11732 | 0057307] 305 | 50 [0 | 2345 | 375 


Momentos Esforços nos apoios 


FREE IT TI TERI am 


Alem, Pilares Arco Pilar Arco | 
M Mp MA Mg | a X Ya 
AB | = 85 pe AR — 3 +2m +23 | —I3-—-15+21— 24 E — £5. 22-— ag | — 4 + 17418 | — 14 + 17 — 18 — 16 + 19 + 20 
1,2 4,354 | 15,155 dás es 4,877 | — ci 
2,4 — | i— | — E5/150 | 14,008 — — 7,191 — 0,076 


(a) (b) 


22576 58 22,576 
Fig. 23 


nula 


m 


a o NI a 
4 
” 
o q a 
- dando codbid. 
» 
” 
. 


IV APLICAÇÃO 


Forma dimensões e solicitação 


| 

3Tn/m.l. 

MODERNO RNDRHA 
B 10 


3STryml. 
NES NanaRIna A 


« 


Fig. 24 


Dada a simetria da solicitação e de estrutura, considerar-se-á apenas metade desta no estudo 
que se segue: 


o em, Secção I 
Pila. Viga. Are. (mº) Eds (m') CPn/m?) (m) 
1,2 mm coma 0,30>x0,30 | 67>< 1074 | 21>< 10 4,00 
3.4 ca | 0,40>0,40 | 21,3 1074 | 21>< 109 4.09 
56 Eu Gi 0,40>0,40 | 21,3x I0-4 | 21 >< 10) 4,00 
7,8 ns sa 0,40 x040 | 21,3xX 1074 | 21>x<I0º 4,00 
9.10 — | — | 0,90x030 | 67x 10-4| 21x 0º 4,00 

aa 2,4 | — 0,30 >x0,50 39L2>< I0"* | 21>< 10º 8,00 
a do | = 0,30><0,50 | 91,2>< I074 21 x 10º 8,00 
asa | -— | 45 0,30 >< 0,40 16 0>x< 1074 21 >< 1G* 12,00 
-— — | 6.8 0.30> 0,40 | 160x 1074 | 21X IO? | 12,00 


Condições de apoio 


Encastramento perfeito de todos os apoios. Os apoios 3 e 7 sofrem assentamentos elásticos. 
Por assentamento elástico entende-se o assentamento que determina uma reacção que lhe é propor- 
cional. Designando por: 


à z—o assentamento (m). 
q—a carga por unidade de superfície que determina o assentimento anterior (Tn/m”) 
R — coeficiente de proporcionalidade (vulgarmente chamado coeficiente de balastro ou 
de assentamento) (Tn/m'). 


Teremos: 
Fazendo À z == Im teremos: 


Cop R e admitiremos R = 2.0C0 Tn/m' 
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o=ci>xg9fz— gx SE —ghithr — t'alt AI c'tig: 


o = + x (g'g£g + g'€gz) + g'Ego'r + ELlgr + g'LeS + S'Eig + g'las: 
o= ox gibr4gxt'alhA-hb>=<ELLgr + riSg + Sgfaz— E'zofor — o'gfor —: 


o=g9xcu'tig + t>g'lzS 4S lol — o'gfgr — o'fgg'b —: 


ogôpIado op Sa10J92] 9 PIDUQN|IUI OP S9JuU9/94309 SOp OBÍLIIIGA 


gte | S'tSi o o o 
[o | [o 9'L6L | g'ntg n'Ega 
'gati |c'trg—| g'£go1 |[ELggi—| O 
o Ss tig— | o | o g'LzS — 
“X | o $ 94 | o! 
EI EI | OL é g 


OgÓBIadO 9P SaloJILA 


(E) 
(=) 


gIFI ic o o Vi Lr's — o b'zlt | ='FIg I co + 
 gEgo | E'LL9T g'LeS o c'o06'1 — | S'EIg g'LeS 1I= a E 
r'1£9 |S gftsc—| O bel S€1g E cobol = | o'gfgI— 1=t1 | E 
o o SColL — | ='hig g'iz o'gtg'1 — |Otggb— |1=& | 1 
E E | = = —— —— — DO 

'N IN IN a bis | IN EN | 
| ' odil | wN 

L | q ç + Fa I | 


| 


11Cg E'ESSz | gvi 
o'g£fg1 | o'glz E chIg | 
S'gfa'z GER | = 
S'£o! o'Lot-1 E 
ENE 
14] | Ás] 


| | | | 


g'tgor | gEz | g'L6L | coz |, orxo'gr| «orx z 
= SEM Ei COM |y OLX GÊ, «OI Iz 
ELEGE | = g'ntg 007  |p-0I>CE'IE| OX Ig 
| 
g'LeS = gtgs | coob |, 01X L'9 | «OLX Iz 
o O a E A 
1X | (uu) | (uu) | UT 
= | [ss | Ia 
8] | ! | I | E 


BIQUBNHuI 9p S9jualdy909 
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Reacções nos extremos dos elementos considerados perfeitamente encastrados 


GC 
E Vie o DOS pe tÔ DOI TA 
2 


vigas 2,4 e 8,10 
É 3x 64 


12 12 


= 16,000 Tnm 


Ve == Vy = 0,500 ql == 0,5>x2><12=-12,000 Tn 


E) 

a >< 
arcos 4,6 e 6,8 Hips E que 2 SS 7 000 Th 
| 8 f 8x2 


Fig. 25 


Acções desequilibradas nos nós 
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Relaxação 


SO E 3 VS en a ante room 
| valor I 2 3 | 4 | 
Tipo da | 
| relaxação No Ny | XA YB 
16,000 — 16,000 — 18 000 12,000 
3 —- 000947 | — 4996 ' — 5810 | 17995 | 0000 | 
| 11,004 — 21,810 | — 0,005 12,000 
2 | — oo021t | 3.456 | 21,738 | — 1,294 | — 0.997 + 
na: 1 4460 | — 0,072 | — 1,299 | — 11,003 
I | 000308 | — 14.424 | — 5.045 1,625 | | 1,892 
RE A ias 0036 | — 5.117 | 0.326 12,895 
4 o c0593 3,642 | 2801 | 0,000 | — 12,875 
di 36798 | — 2,316 0,326 o 020 
I 000078 | — 3653 | — 1,278 | ogtr 0,479 
; |” o025 [— 3.594 | 0737 0,499 
2 — 0 00034 O 55) 3 503 — 0208 | — o161 
| o 582 | — o,ogt 0,529 0,338 
3 O 00027 | O 142 | 0.166 -— 0,513 0,000 
Co | 9) SH) ua " “0,338 
A J  Oo0is4 | — Ojar | — oiase |  cobr | oog4 
| 0,003 | — 0,177 0,007 o 432 
-A |  Socodo]  eiio |  GMbs | | COMO | — 59 
O 113 | — 0,092 O 097 0,041 
I O 000024 — OII2 | — 0039 0,013 0,015 
= sacia 0,001 | — O 131 | oro |  oco56. 
2 — 0,0000112 0.020 | 0,124 | — O 007 — 0,006 
Ra e — 0007 |  o103 0,050 | 
3 0.C000S o 028 0033 | — 0103 O 000 
Co | 0049 | oo26 | 0000 | cos. 
I O 000010 | — 0047 | — ool6 | 0,005 0,006 
0,002 0010 | 0,005 0,056 
4 0,000C07 0,004 0,003 | 0000 | — 0,015 
a 1,006 | 0013 | 0095. OO41 
2 oO 000001 | — 0002 | — o,I0 | O 001 0.000 
E U,004 0003 | 0,006 | oo4 
4 0,000018 O.OLI 0,008 O 000 — 0,039 
(a 0015 | ooftl 0006 0,002 
I 0,0000039 | — 0,014 | — 0,004 0.002 0,002 
Ro 0,001 9,007 | 0008 |  o004 
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Determinação dos momentos e esforços nos extremos das peças 


Elem. [rx] as xx] As (U, — Us) |ry las |yy]4B (Ve — VA) [xrlas [yr] ap E: AB eleclan FA rg | (1.3) (2.3) | (4.6) (5.6) (7.11) (7. 12) (8.11) (8.12) | (9.11) (9.12) (10.11) | (10.12) 
ou al | 
Ab (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (1x) (12) (13) (14) | (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) 
e | e ii) ci iii js | e | io] a | o | pr | per e | see ——— sas "ss E [ES E 
12 527,6 | 263,8|/—0,009146] — — | o 527,60 | — I.407,0 | 703,5 | o 0,004051!| — gba — apt — — o 2,137 | — | — | o 5,709 o 2,850 
34 | 1.677,3|] 838,4 |—0,009146) — — o 1.677,3| — 4.473,0 | 2.236,5 o — 0,002461 — 15,340, — 7,670 — — o — 4,128 -— -— o — II,OIO o — 5,504 
2.4 —— — — 614,2 | 1535 | 0,0061935) — 614,2 | 3.276,0 | 1.638,0| 0004051/—0,002461) — — 3,768 0,942 — -—- 2,488) — r,gt1! 13,271/— 8,062 6,635] —4,031 
4,6 1.063,8 | 797,8 —0,009146| 141,8 23,6 |—-0,006135| 1.063,8 | 141,8 | 2.553,3 | 851,1 |- 0,002461 o 9,729] 7,297] 0,870 |—o,145| — 2,618 o — 0,349 o — 6,284 o — 2,094 o 
| 
Momentos Esforços nos apoios 
Elem. Pilares Vigas | Arco Pilares no Vigas Arco 
MA M, MA My MA Mg A EA E; | É 
AB |-I3+-2e +24 |—-134-22+23|)-15+214+24 | —15+224+23|— 13 +15+21—24 |I3—-15+22—23| —14 + 17+418| — 16419-420 — 14 +17—18| — 16 + 19420 
enmeuaos,-zzzmtmá8õãõãíãÇí888 E EEE ——————, | e À eee | caem 
ta 7,675 10,525 — -— — — | 4,549 Eis E RE 
3,4 9,836 4,330 eme e — — | 3,542 — Eee E 
24 = -— 5,472 — 5195 — — — 0,035 — [0 — 
te 
4,0 poi E E | pao — 15,143 12,653 Ea ci — 9,915 — 0,204 
12,000 12,000 12,000 12,000 


10525 


4594 


12,035 
6805 108 


| 


11,965 24,239 


POSTOS DE TRANSFORMAÇÃO 
BLINDADOS 


PARA INTERIOR E EXTERIOR 


Posto de transformação de 250 KVA, 6,6 KV, para interior, 
fornecido para os Serviços Autónomos de Luz e Água de Luanda (LAL) 


CONSTRUÍDOS POR: 


JAYME DA COSTA, LDA 


LISBOA — PORTO — LUANDA. 
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A humidade éa causa de sérios 
problemas de ordem técnica e de 
incalculáveis prejuízos. 


Estes inconvenientes podem ser 
totalmente eliminados com a uti- 
lização dos Exsicadores 
“BIALEC” fabricados por 


BIRLEC LIMITED 


BIRMINGHAM — Inglaterra 


Representantes para Portugal: 


Ohlers, Lindley, PLinitada 


Rua Ferregial de Baixo, 33-2.º Palácio do Atlântico, 408 
LISBOA Praça D, João 1 
Telef. 2 1321/4 PORTO — Telef. 31414/5 
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e e 


Verificação da relaxação 


Valor I 2 | 
da E ra N, N, | 
16,000 | — 16,000 | — 18,000 12,000 
I 0004051 | — 18,971 | — 6,634 |— 2137] 2,488 
2 | —oooaór | 403t| 25355 — 1508 | — 1160 
3 — 0,009146 | — 4,826 — 56 | 17.379 | 0,000 
RE. 1» 0,006135 | | 3.767. 2,894 E 0,000 | — 13,320. 
EA 0,001 0,007 | 0,008 | 0,004 | 


Valores totais das rotações e translacções 


Ta == 0,004051 rad. 

r, == 0,002461 rad. 
U, = 0,009146 m 
V; = 0,006135 m 


Verificação das reacções de apoio 
Dos valores obtidos na fig. 27 resulta 


4,594 — 4,594 + 3,542 — 3,542 =0 


Equilíbrio de forças | 48 + 48 — 11,965 — 22,239 — 23,592 — 22,239 — 11,965 = 0 


— (24 + 36) — (48>< 20) — (24><4) + 11,965 x 40 + 24,239 x 


Equilíbrio de momentos | > 32 + 23,592x 204 24,239x 8=0 


Deverá ainda proceder-se à verificação do equilíbrio de forças e de momentos nos nós. 


NOTAS INFORMATIVAS 


C.D. U. 624.344.5/94 


Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 


— Elementos extraídos das estatísticas mensais do Repartidor Nacional de Cargas (R. N.C.) 


Nora: 


As produções e os consumos das empresas do R, N. C. representam 


cerca de 91,6 */, dos totais do Pais, 


SETEMBRO 


I — Breve nota mensal 


Do ponto de vista hidrológico o mês de Setembro foi 
um pouco inferior ao normal. 


II — Elementos gerais (GWh) 


a) Mensais 


) 1956 | 1957 Eds 
0 
Produção hidráulica (Ph). 149,6! 1451 | — 3 
Produção térmica (Pr)... ... 0, p| (1,0 (Ú 
Produção total (PT)... ... 149,6] 1451 |— 3 
Cons. electroquímico (Ceg) (1) | 30,4 1,1 |— 63 
Cons, permanentes (U + (1) | 112,0 1260 +. 13 
Consumo total (CT) . ... (1) | 1424| 1371 ER po 4 


b) Acumulados desde 1 de Janeiro de 1957 


Variação 
1956 | 1957 | oyo 
u 
Produção hidráulica (Ph). 1446,8! 14586 + 1 
Produção térmica (P+), . «+. 36 51,1 - 
Produção total (PT). - o .| 1450,4 1510,2 + 4 
Cons. electroquímico (Ce (1) | 350,1) 266,6] — 24 
Cons. permanentes (Cp). (| 1048,7 | 1167,9 |+ 11 
Consumo total (CT) . +... vd Pa ri 1395,8 | 1434,5 | + 3 
Nota : 


(1) Vidé nota referente ao mês de Janeiro de 1957, 


II — Diagramas de carga dos dias característicos 


4.º feira; 
19-9-956 | 18 9-957 
Produção hidráulica (Ph) — MWh| 5518 5268 
Produção térmica ( P+)— MWh. . O 0 
Produção total (Pr) — MWh ...| 5518 | 5268 
Utilização da ponta (U) — horas 17,2 | 120 
Factor de carga (a) . . «+. 012 | 0,71 
E 
Mid o. db) db 
Pot. máx. 
TÉCNICA 
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IV — Energia armazenada nas principais albufeiras 
no fim do mês. 


Energia armazenada 
Albufeira Po EM 


GWh | 0% (1) 
Venda Nova « vc « «vas 29,1 23,2 
Salamondo . . «cv cos. 19,3 69,9 
CANIÇÕCA «cocos wc. 8,6 25,9 
Quilhofral e sá mas sc 2,2 26,5 
Lagoa Comprida + . .... 8, 29,0 
Santa Luzia « «cw sw ds 22,4 96,4 
CABEM dm snes us 64,1 18,9 
Castelo do Bode. . .... TT, 47,5 . 
Pracana cocos. 1,5 121 
4,6 (2) 26,9 

238,6 28,8 


POVOÃ o uz E vm a j E tá 


Notas: 


(1) Coeficiente de enchimento em energia das albufeiras, 
definido pela relação 


Energia armazenada 
a memo DOT Og 
Máx. energia armazenável 


(2) Inclui 1,6 GWh armazenados no açude do Poio, 


TEODOLITO TAQUEOMÉTRICO FTI 


Ota asa 


MODERNO APARELHO DISPONDO DE: 


O Observação simultânea dos 
dois limbos pelo microscópio 
colocado so lado da ocular da 
objecliva 


O leitura por sistema directo de 
lc e décimos por eslimaliva 


€ leitura do nível vertical por 
coincidência 


€ Prumo óptico 


€ Grande luminosidade 


Características ópticas e mecânicas 


- 


Telescópio analóciico com lente de focagem Nível junto so círculo vertical legível por ajusta- 

interna, retículo medidor de distância, telescópio mento da coincidência 40". Passo circunferen- 

reversível com o extremo da objecliva. cial sobre aa de “io 360º ou 4009, 
elooudelo. 


Conslante de adição. . . 
Constante de multiplicação . 
Comprimento do telescópio. 
Aberlura da objeclivo +... 
AMPRAÇÕO cs s.o o So 0x4 
Dislôncia mínima de mira . . 
Nível redondo. «. «cv... 
Nível de reversão. . «+. 
Nível de horizontalidade paralel 
eixo de Inclinação. +. « . «cv. 


100 Diôêmetro do círculo horizontel 90 mm, do clr- 
190 mm culo vertical 70 mm 


36 mm leitura directa horizontel 1º fe 
28x Leitura por eveliação « . 1º le 
Mo Leitura directa vertical . 0,1 = 6º 0,2 e =920cc 
6a 8 Leitura por avaliação . . 0,2 == 12" 0,25€ = 25cc 


. 30" 
o Pesos: Instrumento 4,7 kg, Tripé com pés exten- 
40” síveis 6,0 kg, Caixa metálica 4,0 kg 


. 
(6) 


REPRESENTANTES 


PAPELARIA FERNANDES 
L. DO RATO, 13—TELEF. 6611167 LINHAS —R. DO OURO, 145— TELEF. 28361 


LISBOA 
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Nas Indústrias do Petróleo 


Desde a extracção 


O! I 
dy 


ú ii Tm: 
pu e 


é conhecido em todo o mundo 


Distribuidores exclusivos: CANELAS & FIGUEIREDO, L.ºA 


R. dos Fanqueiros, 46 —- LISBOA 
2 5058 
Telef. 4 24502 
2 1729 Teleg. CANELASCO 
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C. D. U. 550.834 


PROSPECÇÃO GRAVIMÉTRICA 


O Problema da Interpretação 
PELO ENG. DE MINAS (1.S.7.) JOSÉ QUINTINO ROGADO 


(Continuação) 
CAPÍTULO II 


A indeterminação da configuração 
potenciante 


Temo-nos até aqui referido apenas ao campo 
gravítico, que pode pôr o prospector em con- 
tacto com as ocorrências geológicas, dando delas 
uma imagem, manifestação física da sua massa. 
Ocorre agora perguntar se será possível, par- 
tindo do conhecimento desse campo, estruturar 
univocamente as configurações geológicas por 
ele responsáveis. 

Infelizmente a resposta é, como se sabe, ne- 
gativa. Não é possível estabelecer tal correspon- 
dência unívoca. Quere dizer, a interpretação dos 
dados gravíticos em termos de massas poten- 
ciantes esbarra com uma determinação intrin- 
seca, o que tem dado lugar ao desenvolvimento 
do que se pode chamar «processos de interpre- 
tação condicionada», os mais elaborados dos 
quais partem precisamente da análise daquela 
indeterminação. Torna-se por isso conveniente 
que anteriormente ao estudo da interpretação se 
faça em termos precisos a definição da indeter- 
minação. 

A análise de Green à função potencial per- 
mite, como apontamos, considerar o potencial 
gravítico U (P), na superfície da Terra, devido 
a qualquer distribuição densitária subterrânea, 
como sendo o potencial provocado por uma dis- 
tribuição de simples camada de densidade 


dg (pm o | ) 
/ Q 


“q 


2xk 27% k ez 


existente sobre um plano horizontal arbitrário S' 
de ponto corrente Q, situado entre a superfície 
da Terra e a distribuição potenciante real. Isto 
é, como se escreveu em (41-1): 


DEM (a 
U(P)= 5 7. : g (O) dS" (41.1) 


Bolseiro do 1. A. €. 


Denominaremos qualquer distribuição subter- 
rânea arbitrária de simples camada, satisfazendo 
à condição de ser 


de «camada equivalente». 

Pois bem, como, uma vez determinada a fun- 
ção potencial, tudo fica conhecido no campo 
dela derivado, segue-se que se removermos a dis- 
tribuição potenciante real e a substituirmos por 
uma camada equivalente nada se alterará naquele, 

Ao plano S”, suporte da camada equivalente, 
apenas se exigem as condições de ser paralelo à 
superfície da Terra e de estar situado entre 
aquela e a distribuição potenciante real. Tra- 
ta-se, efectivamente, de condições muito fracas, 
que se traduzem na existência de uma infinidade 
continua (1) de camadas equivalentes. 

Ora é precisamente nisto que reside a inde- 
terminação da configuração potenciante. 

Do ponto de vista da gravidade essa indeter- 
minação traduz-se pela seguinte equação que 
resulta de (41-I) por simples derivação se- 
gundo a vertical, feita no ponto P 


JU 1 d /1 á 
8 (= ( j io 1 Ea) (O) do 
(1-1) 


e que é uma equação integral pois g (Q) não 
é conhecido. 

A integração de (1-Il) teria a vantagem de 
relacionar a densidade da camada equivalente 
com os valores da gravidade observados, permi- 
tindo construir aquela, o que, além de tornar bem 
patente a indeterminação, forneceria um pro- 


(1) Entenda-se com a potência de um conjunto conti- 
nuo. Como se sabe, a análise de Green parte até de con- 
dições mais gerais. 
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cesso de interpretação, desde que se conheça a 
profundidade da estrutura potenciante real. 


Assim, sendo U o potencial, será Lap U=O, 
cuja solução se pode pôr na forma [25]: 


p—— 
vs sen mxsenny —VYmitn?,z 
E cos mx cosny 


(2-1) 


onde z é a cota em relação a um plano de refe- 
rência, por exemplo, a camada equivalente. 
Nesta, como já vimos, será: 


doU 


Pois bem, se admitirmos que a densidade dis- 
tribuída sobre a camada equivalente é expremí- 
vel por uma função que pertença à classe das 
funções com desenvolvimento em série Fourier 
poderá escrever-se: 


(3-1) 


sen mx sen ny 
cos mx cos ny 


ed (x, y) = 2 à Ban (4-I]) 


e, portanto, por (3-1), será: 


Mas g (x,y, z) é também harmónica e por- 
tanto a solução de Lap g =0, no plano da super- 
fície da Terra suposto à cota + h, poderá por-se 
na forma: 


sen mx sen ny 
mn cos mx con ny (5-1) 


g (xy +h)=22C 


Mas, sendo 


<a 
= t , h 
- angu g (x, y, h) 


/ 


sera, por (2-1) e (5-Il): 


ou 


que é a relação procurada e mostra ser a densi- 
dade da camada equivalente função da profun- 
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didade h a que se encontra, o que é mais uma 
maneira de traduzir a indeterminação que vimos 
referindo. Convém notar que em (7-1I) se supoz 
que a densidade da camada equivalente pertencia 
a uma determinada classe de funções e, portanto, 
não engloba todos os casos possíveis, embora 
a classe daquelas funções seja muito extensa [26]. 

Pode, contudo, argumentar-se que camadas de 
densidade variável não são geológicamente pre- 
sumiveis e, portanto, não só não há interesse 
prático em referi-las, como também o facto de 
poderem substituir qualquer massa potenciante 
não representa embaraços para quem interpretar 
os dados gravimétricos, visto que as não consi- 
derará [27]. 

Ora, notemos que a densidade superficial 
variável pode traduzir ondulações de uma 
estrutura, pois se puzermos 5 (x, y) = y At, da 
onde y é contraste densitário e dt a espessura 
camada e, se admitirmos y é constante, a espes- 
sura será variável — Ààt (x, y). Poderemos pois 
obter estruturas equivalentes como as da Fig. 6. 
Isto é, estruturas com o mesmo campo que a 
configuração real e já presumíveis sob o aspecto 
geológico. 

Além disso, sempre que a camada equivalente 
apresentar duas zonas contíguas com densidades 
de sinal contrário, poderemos transformá-la numa 
estrutura com um único contraste densitário como 
o da Fig. 7-a. Basta considerar que uma tal estru- 
tura, se for extensa e limitada por um plano 
horizantal não produzirá anomalia gravítica local. 
Esta, se existir, resultará apenas da ondulação 
da superfície que determina o contraste den- 
sitário (V. Apêndice) e, nesse caso, o campo 
anômalo será igual ao produzido por uma camada 
equivalente com uma área de densidade pcsitiva 
contígua a outra de densidade negativa, visto 
que, quanto à anomalia que produz, aquela 
estrutura se pode assemelhar a uma estrutura 
como a da Fig. 7-b [25], onde a altura z acima da 
camada equivalente se pode aproximadamente cal- 
cular [27], desde que a estrutura seja suficien- 
temente profunda, pondo 


g (0) | 
27k.y é 


onde y é o contraste densitário admitido. 
E, mesmo que a camada equivalente apresente 
só uma zona de densidade positiva, ainda é pos- 


sivel substituí-la por uma estrutura equivalente 
com um único contraste densitário, por um 
processo análogo ao anterior (Fig. 8). 

Vê-se assim que a noção de camada equivalente, 
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sificadas de regionais e em oposição, no segundo 
caso, chamadas locais. 

O prospector poderá pretender conhecer qual- 
quer desses dois grupos de ocorrências e conforme 


V/liit!i) camada equivalente 

GAS 

eme esfrulura equivalente 
Fig. 6 


extraída da análise de Green, embora traduza 
a indeterminação da configuração potenciante, 
permite construir hipóteses sobre a mesma de 
grande valor para técnicas interpretativas [25], 
o que, se constitui vantagem evidente, é demons- 
tração perfeita da importância da análise da 
indeterminação naquelas técnicas. 


CAPÍTULO III 


Anomalia regional e anomalia local 


Feita a correcção de Bouguer e a correcção isos- 
tática, se necessário, a anomalia de Bouguer tra- 
duz, como já referimos, as anomalias na distri- 
buição densitária subterrânea, isto é, as ocorrên- 
cias geológicas cujas densidades se não coadunam 
com as hipóteses feitas nas correcções. 

Estas ocorrências geológicas podem ser pro- 
fundas e extensas ou, pelo contrário, ocorrerem 
junto à superfície da Terra e terem dimensões 


reduzidas. No primeiro caso, costumam ser clas- 
aquele sobre que incide o seu interesse assim 
estará realizando prospecção regional ou local, 
respectivamente. Compreende-se que, por exem- 
plo, ao realizar prospecção local, haja que se- 
parar a amomalia devida às ocorrências regio- 
nais, que não interessam, das outras provocadas 
pelas configurações locais que se procuram. 
E inversamente, na prospecção regional. 

Surgem pois estes dois conceitos difíceis de 
definir, aos quais já fizemos referência no Capí- 
tulo I, e que traduzem uma necessidade prática 
ineludível — a anomalia local e a anomalia regional. 

Este aspecto da apreciação das anomalias, 
posto desde os primeiros tempos da prospecção 
geofísica, consiste em definir duas componentes 
do campo e atribuir-lhes causas geológicas de 
tipo diferente, problema que, se atendermos ao 
que ficou dito sobre a indeterminação da massa 
potenciante, é insolúvel, no caso geral. 

É claro que o conhecimento geológico regional 
ou local da área a interpretar é de incomparável 
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valor, porque pode conduzir a situarmo-nos 
apenas na consideração de alguns casos parti- 
culares mais prováveis. Mas, infelizmente, os 
dados geológicos iniciais são cada vez menos 


direcções à custa de qualquer conhecimento geo- 
lógico estrutural ou por simples verificação de 
determidada lei que ressalte do conjunto dos 
dados geofísicos de que se dispõe, nos quais 


superfície da Terra 


contraste 


7 densifário 
7 posilivo 
PAPA AP 
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(a)-Eslrulura com um único contraste densilário 


Superfície da Terra 


contraste 


4 


densilario posilivo 
Ei densidade posiliva 


densidade megaliva 


DMI III III = —= 


contraste densi- 


(b)-Eslrulura equivalente sob o ponfo de visfa da anoma- 
lia de Bouguer e camada equivalente 


Fig. 7 


frequentes ou mais escassos, particularmente os 
dados geológicos locais, e por isso a prospecção 
geofísica tem, as mais das vezes, de lançar mão 
de métodos próprios que lhe permitam delinear 
o que é regional e o que é local, nas anomalias 
observadas. 

Na essência todos esses métodos se resumem 
em realizar uma atenuação da anomalia observada, 
entendendo-se por tal a operação que suprime 
certos aspectos da anomalia, fazendo sobressair 
outros. 

Assim, a atenuação pode ser conduzida de tal 
modo que os valores altos do gradiante hori- 
zontal da anomalia, segundo uma ou mais direc- 
ções, sejam suprimidos, tendo escolhido aquelas 
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certas configurações da anomalia terão então 
uma «permanência» que permitirá deduzir essa 
lei, enquanto outras serão excepcionais. Far-se-á. 
portanto, com atenuação deste tipo, a ablação 
das anomalias locais e ficar-se-á de posse da ano- 
malia regional, à qual se associou a ideia de ser 
provocada por configurações geológicas profun- 
das ou extensas, o que está, na verdade, de 
acordo com o conceito habitual de regional. 
Noutro tipo de atenuação, procura-se suprimir 
as anomalias locais impondo uma representação 
analítica simples (um polinómio inteiro de baixo 
grau, quase sempre) à anomalia regional e deter- 
minando os parâmetros daquela representação 
por processos de interpolação, a um dos quais 


já nos referimos no Capítulo I. Como se vê, não 
se trata de forma de atenuação muito diferente 
da anteriormente referida, porquanto escolher 
um polinómio de baixo grau para representar a 
anomalia regional é impor determinados valores 
aos gradiantes horizontais daquela. 

Ainda um terceiro tipo de atenuação, que na 
cronologia não é o último, procura-se suprimir 
a anomalia local, considerando-a análoga a uma 
dispersão sujeita a qualquer lei de probabilidade. 

Evidentemente, o princípio que se adopta, 
para construir a anomalia regional caracterizará 


uma área extensa, e pessoas diferentes obterão 
os mesmos resultados, ao realizarem a atenua- 
ção. Mas, se a estrutura geológica regional é 
complexa ou não é bem distinta das ocorrências 
locais, o resultado da atenuação será seriamente 
afectado pelo método escolhido e, em conse- 
quência, por elementos pessoais de quem fez a 
escolha, 

Da verificação destes factos resultou o apare- 
cimento de numerosos métodos para separar a 
anomalia local da anomalia regional. Podemos 
reuni-los em dois grandes grupos — métodos grá- 
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Fig. 8 


o método a empregar e os resultados que se 
alcançarão. E não é demais frizar a importância 
daquele na interpretação dos dados geofísicos. 
Por isso, esta escolha da forma como se irá 
definir as anomalias regional e local deve ser 
criteriosamente pesada e encarada sempre num 
conjunto tão extenso quanto possível de dados 
geofísicos, geológicos, mineiros, etc. 

Em certos casos a separação do que é regio- 
nal na anomalia é imediato, como, por exemplo, 
no caso já citado, da existência de um primeiro 
gradiante horizontal da anomalia mais ou menos 
“uniforme ou variando lentamente ao longo de 


ficos e métodos analíticos. Quase todos eles le- 
vam à construção de dois campos, o campo 
regional e o campo residual, diferença daquele 
para o campo original, e que é considerado o 
campo local, 


1 — Métodos gráficos de atenuação 


Estes métodos apresentam duas modalidades 
que conduzem sensivelmente aos mesmos resul- 
tados, quando forem bem aplicados e os valores 
originais de que partirem forem em abundância. 
A primeira modalidade, a que chamaremos «Mé- 
todo de Atenuação dos Contornos», é a forma 


TÉCNICA 
109 


mais directa de realizar a atenuação e consiste 
em obter o campo regional por regularização 
dos contornos da carta do campo original, en- 
tendendo-se por tal a substituição dos contornos 
mais ou menos curvilineos por outros rectilineos 
ou de pequena curvatura, o que corresponde, 
evidentemente, a admitir que a superfície analí- 
tica que representa o campo regional é bem 
interpolada por um polinómio de baixo grau. 

Este método é muito empregado em razão da 
sua simplicidade, e é suficientemente seguro, 
quando ocorrem configuracões regionais nítidas 
ou o andamento dos contornos residuais for 
diferente do dos contornos regionais. 

A outra modalidade que denominaremos «Me- 
todo de Atenuação dos Perfis» é já mais evo- 
luída, mas também mais trabalhosa, e resume-se, 
na sua forma mais simples, em construir perfis 
do campo original, paralelos e equidistantes, e 
cuja direcção comum ou é mais ou menos per- 
pendicular aos contornos do mapa original ou 
paralela a qualquer direcção geológica previle- 
giada, perfis esses que são depois regularizados 
suprimindo as zonas de curvatura acentuada e 
procurando manter, tanto quanto possível, o 
mesmo gradiante horizontal ao longo de cada 
um deles. Finalmente, à custa dos novos perfis 
assim obtidos, constroi-se o campo regional e, 
por diferença para o original, o campo residual. 

Numa forma mais complexa, este método rea- 
liza-se segundo duas direcções, ortogonais ou 
não, de modo que se obtém, em cada ponto do 
cruzamento de dois perfis, dois valores que 
devem ser muito próximos para se tomar a sua 
média como o valor do campo regional nesse 
ponto. Se assim não for, há que repetir em 
ambos os perfis a atenuação até os dois valores 
obtidos em cada nó diferirem pouco. É uma mo- 
dalidade aconselhável quando o campo original 
se apresenta complexo. 

Todos estes processos gráficos, para serem se- 
guros, têm de ser efectuados por pessoa mais 
ou menos conhecedora do significado dos concei- 
tos em que se estruturam. Apesar disso, ou talvez 
por isso mesmo, os resultados a que conduzem po- 
dem ser muito influenciados por factores pessoais. 


2 — Métodos analíticos de atenuação 


Com a intenção de automatizar a atenuação e 
assim eliminar o principal defeito dos métodos 
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gráficos, estruturaram-se vários métodos anali- 
ticos que se podem classificar em «Métodos de 
Atenuação pela Média» e «Métodos de Atenua- 
ção por Interpolação», embora todos eles sejam 
baseados em operações de interpolação [32]. 
Qualquer deles, por não exigirem especialistas 
para a sua execução, embaratecem as operações 
de atenuação e são, por essa circunstância, bas- 
tante usados em trabalhos de rotina, 


2.1 — Métodos de atenuação pela média 


O mais simples método de atenuação pela 
média consiste em calcular as médias aritméticas 
dos valores da anomalia original que caem den- 
tro de uma área, quase sempre quadrada, que se 
centra sucessivamente nos pontos, onde se pre- 
tende conhecer o campo regional, e considerar 
em cada um desses pontos este campo como 
sendo igual às médias obtidas. 

É costume calcular por interpolação o campo 
original em pontos regularmente localizados na 
área prospectada, antes de se efectuar esta ope- 
ração. 

Como é evidente, o valor do campo regional 
é muito afectado pela área escolhida para o cál- 
culo da média. Se esta é extensa, o regional ten- 
derá para valor baixo e o campo residual conse- 
quentemente aproximar-se-á do campo original, 
Se é pequena, o campo regional aproximar-se-á 
do original e o campo residual tenderá para 
ZeTto., 

Outra forma simplista de atenuação analítica 
consiste em considerar sobre a área prospectada 
uma rede de malha regular, definida por duas 
direcções ortogonais, calcular por interpolação o 
campo original nos seus nós e construir, ao 
longo de cada fiada da rede, as médias aritmé- 
ticas de um certo número de valores do campo 
original simeétricamente dispostos para um e 
outro lado de cada nó. 

A cada nó da rede ficarão, pois, correspon- 
dendo dois valores médios, um por cada direc- 
ção da rede. Toma-se por valor do campo regio- 
nal a média aritmética desses dois valores. 

Como é fácil de mostrar, este modo de proce- 
der pode considerar-se baseado na hipótese de o 
campo regional ser interpolado por um polinó- 
mio inteiro de qualquer grau, e não apenas do 
primeiro, como tentou fazê-lo Agocs [32]. 

Além disso, porque se considera para o cál- 


culo do valor regional num ponto os valores do 
campo original apenas segundo duas direcções, 
este método pode conduzir a resultados erróneos. 

Quer este último método, quer o anterior, 
podem considerar-se processos empíricos de reali- 
zar a atenuação. Serão, quando muito, aconse- 
lháveis em casos favoráveis. Consegue-se com 
eles automatizar a atenuação mas obtêm-se resul- 
tados nitidamente mais grosseiros que com os 
métodos gráficos. 


2.1.1 — Método de atenuação com grade 


Sem dúvida, o processo de atenuação que tem 
sido mais usado, quando o volume de dados a 
estudar é já elevado, é aquele que dominaremos 
«Método de atenuação com grade». 

Com este método pode obter-se directamente, 
quer o campo regional, quer o campo residual. 
Interessa-nos especialmente abordar o segundo 
caso. 

A forma mais simples de construir o campo 
residual em determinado ponto, utilizando aquele 
método, consiste em calcular a média dos valo- 
res do campo original situados sobre uma cir- 
cunferência centrada nesse ponto e deduzi-la ao 
valor do campo no centro. Covém evidentemente 
que os pontos tomados na circunferência, para 
construir a média, sejam equidistantes e no maior 
número possível. Por isso, se os pontos obser- 
vados se não distribuem regularmente, há que 
fazer essa distribuição por interpolação. 

Nas aplicações usam-se para o cálculo das 
médias 4, 6, 8 ou mais pontos. Para facilitar a 
escolha dos pontos desenha-se em papel trans- 
parente (1) e na escala da carta da anomalia de 
Bouguer, uma grade, onde estão marcados o cen- 
tro e os pontos onde devem ser feitas as leituras. 
Daí o nome deste método. 

O campo residual construído com grade de 
uma só circunferência depende considerável- 
mente do diâmetro daquela e do número de 
pontos tomados para o cálculo da média. 

Por essa razão se procurou modificar o método, 
utilizando, em vez de uma, várias circunferências, 
construindo a média ponderada das médias 
aritméticas obtidas em cada uma delas e dedu- 
zindo-a ao valor central para se obter o campo 
residual. 


(1) Ver Capítulo 1, Fig. 5. 


O modo de estabelecer os pesos dessas médias 
é arbitrário, embora se possa procurar dar-lhe 
uma estruturação analítica, que elucidará sobre 
o significado do sistema de pesos empregado. | 

Antes de mais, note-se que, usando uma única 


circunferência, o campo residual À g é definido 
do seguinte modo: 


à g=g (0)-—g (r) (1-1) 
onde g (0) é o campo original na estação e g (r) 
a média na circunferência de raio r. Portanto, 
como vimos, será: 


g()=ao+tartart+... 


Supondo, como já fizemos no Capítulo I, que 
g (r) numa pequena vizinhança é bem aproxi- 


mado por 
G(r)=açãas rº 


teremos : 
g=g (0)-aç—as rº 
Ora: 
ão = g (0) 
donde: 
= — 24 Tº (2-1) 


Mas, já vimos (Capítulo I) que, se az tiver 
sido calculado numa vizinhança da origem eva- 
nescente, será: 


E por isso pode escrever-se: 


-p? q? 
Fa g == + (—s é 
4 d z* r=0 


Quere dizer, o valor residual é, a menos um 
factor, o valor da segunda derivada, qualquer 
que seja o raio r onde se calculou a média, 
desde que o valor de az escolhido esteja cor- 
recto numa vizinhança da origem arbitráriamente 
pequena. É claro que, para as ser válido numa 
vizinhança evanescente e numa vizinhança fi- 
nita, há que admitir que a aproximação G(r) de 
g(r) não é uma aproximação, mas sim o valor 
correcto, isto é, G(r) = g(r). 

Portanto, introduzindo a hipótese de g(r) ser 
um polinómio do segundo grau, o campo resi- 
dual construído do modo referido é, a menos 
um factor, o campo do segundo gradiante ver- 
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tical da anomalia. Mas para que g(r) seja um 
polimónio do segundo grau, terá que ser a ano- 
malia um polimónio do segundo ou do terceiro 
graus em r. É esta, pois, a condição para a 
identificação do segundo gradiante vertical com 
o campo residual, qualquer que seja o raio r de 
circunferência onde se calcula a média. 

Quando aquela condição não for satisfeita, o 
campo residual não coincide com o campo do 
segundo gradiante. Chama-se-lhe por isso campo 
derivado [33]. 

No caso geral será, portanto: 


e(t)=a+ar-tar'+o... 


e se adoptarmos, para definição do campo resi- 
dual, a expressão (1-III), importa calcular sem- 
pre g(r). 

Já vimos no Capítulo 1, como calcular gráfica 
ou analiticamente o coeficiente az. O cálculo 
de aí, ... far-se-ia de forma análoga; contudo, 
não se costuma fazer. 

Quere dizer, habitualmente calcula-se o campo 
residual, usando o processo adoptado para o 
cálculo do segundo gradiante vertical. Este modo 
de proceder estrutura-se no seguinte [33]: 

A anomalia satisfaz a equação de Laplace 


E, se representarmos os valores do campo, 
nos pontos tomados na circunferência de raio r, 
por b; e o valor central por €C, será, tomando 
quatro pontos situados em dois eixos hortogo- 
nais: 


Pg 1[(bi—C (C-bi)] (bi+b:—2C) 
J x? T r r pe rº? 
Pg 1[(bi—C (C—-b)] (br+bi—20) 
d yº r r r aee r2 = 
e portanto 
Pg 4. em (bi + bz + bs+bi) 
dz tl 4 


Deslocando os eixos n vezes, teremos: 


o? g 4 
dz op 
C— (bi+ba+ba+ bg) +..+ (ban-3+Dan-2+Dbyn-r+ba ) 
4n o! 
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e sen for suficientemente elevado, o parêntesis 
recto do segundo membro representará a verda- 
deira média sobre a circunferência de raio r, 
isto é g(r) e, portanto, o campo residual, no 
caso geral, construído à custa daquela média, 
será também igual (a menos um factor) ao 
campo do segundo gradiante vertical se o raio r 
for suficientemente pequeno. 

Pode, por isso, substituir-se sempre o campo 
residual pelo campo da segunda derivada, esco- 
lhendo r arbitráriamente pequeno, e fazer, por- 
tanto, o cálculo daquele, à custa de as. Isto é, 
mesmo no caso geral, subsistirá (2-IIl) e conti- 
nuar-se-á a chamar ao campo residual, assim 
obtido, campo derivado, quando, por razões eco- 
nómicas, não houver interesse em escolher r 
suficientemente pequeno. 

Ora, como o cálculo de as, como vimos, se 
faz considerando vários valores de r (isto é, 
várias circunferências concêntricas) e interpo- 
lando, para se obter o melhor valor, resulta do 
que se disse, a mesma técnica para o cálculo do 
campo derivado e simultâneamente um sistema 
de pesos para as diversas médias, que vai buscar 
a sua razão de ser às semelhanças entre o campo 
derivado e o campo do segundo gradiante vertical. 

Existem, pois, para expressões de a», além 
das já indicadas no Capítulo I, várias outras 
baseadas na utilização de três ou mais circunfe- 
rências e que podem ser vistas em [21, 24 e 31]. 

Em conclusão, a atenuação com grade estru- 
tura-se na circunstância de que o campo residual 
construído (campo derivado) se aproxima do 
campo do segundo gradiante vertical, o que vale 
o mesmo que dizer que se considera o campo 
regional como um polinómio de segundo grau, 
porquanto a aplicação do operador 


o Ps A 

2º dx? dy? ) 

o suprime. É claro que aquele operador, no 
caso geral, é aplicado de forma aproximada, 
visto não haver conveniência em usar no cál- 
culo de az circunferências de raio suficiente- 
mente pequeno e portanto a sua função de atenua- 
ção [33] é dependente daquele raio. 

Esta definição de anomalia regional, implícita 
naqueles métodos de atenuação, é plenamente 
justificada se atendermos a que se chamou ano- 
malia regional àquela que é provocada por con- 
figurações geológicas extensas ou profundas, e 
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constitui o ponto de partida para os métodos 
de atenuação por interpoloção com polinómios 
inteiros. 


2.2 — Métodos de atenuação por interpolação com 
polimónios 


Estes métodos recorrem à interpolação pelos 
minimos quadrados e definem o campo regional 
como sendo um campo bem interpolado por 
polinómio inteiro de baixo grau. Em oposição, 
admitem que as anomalias locais só podem ser 
bem interpoladas por polinómios do grau elevado. 

Contudo, para a estruturação de um e de outro 
campo falta, nestes métodos, precisar, em cada 
caso concreto, quais os graus dos dois polinómios 
interpoladores. 

Ultimamente [34] surgiu uma tentativa para 
realizar essa definição, que referiremos à frente, 

Em conjunto, este modo de atenuação tem 
cada vez maior emprego em resultado do uso 
sempre crescente das máquinas calculadoras- 
Trata-se, efectivamente, a nosso ver, da maneira 
mais criteriosa de realizar a separação dos dois 
campos, sem deixar de ter a elasticidade neces- 
sária, para se poder adaptar a casos muito 
diversos e a hipóteses muito diferentes. Contudo, 
só se aconselha, quando o volume dos trabalhos 
justificar a estandardização que o seu emprego 
implica. 

Na sua forma inicial [35], este método de 
atenuação punha o problema nestes termos: 

É conhecido o campo original g (x, y) num 
conjunto discreto de pontos dispostos nos nós 
de uma rede de malha constante, que represen- 
taremos por S, situada no plano Oxy coincidente 
com a superfície da Terra. 

Pretende-se interpolá-lo por um polinómio 
inteiro em rey de grau n, fixado, e que será 
a representação analítica do campo regional. 
O campo ou anomalia local é, em cada ponto, o 
desvio entre a interpolação e o valor do campo 
original. 

Isto é, sendo R (x, y) o polimónio interpolador, 
representação analítica do campo regional, será: 


g (x y)=R(x,y+ Ag (x,y) 


onde, como sempre, & g (x,y) é a anomalia local. 
É pois 


Rixy)= 53 


K, |=0 


Cir xK y! 


Ponha-se 
nº (Cx) =S[g (x,y) —R (x, y)]º 


onde representamos por S a soma efectuada com 
todos os valores da rede 5. 

As condições de mínimo, inerentes ao método, 
que permitrrão determinar Cy; serão dadas pelo 
sistema das equações normais 


nn-—K 
3 ESTE yri= Se NAN Y 


, j=0 


A solução deste sistema pode simplificar-se e 
para tanto dê-se-lhe a forma matricial: 


Is srt yi+) k (Cii] E Is g (x, y) xt yih 
cuja solução é, evidentemente: 
(Cu) = (Sxitk yi+j-, IS g (x, y) xt y') 


Ora, se atendermos a que os valores de g (x, y) 
são conhecidos nos nós de uma rede rectangular 
de malha constante, como a da Fig. 9, os elemen- 
tos da matriz 


is xi+k yi+t) 


podem escrever-se : 
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donde se deduz que os elementos da matriz só 
não são nulos, quando (i+k) e (j +41) forem 
ambos pares. Conclusão esta que é geral. Quere 
dizer, utilizando sempre a mesma malha e poli- 
nómios do mesmo grau, para interpolar o campo 
regional, o cálculo do coeficiente C, é muito 
simplificado, embora mesmo assim possa ser 
muito trabalhoso. 

Como se viu, houve que arbitrar o grau do 
polinómio interpolador do campo regional, o que 
pode fazer-se por estudo da carta da anomalia 
de Bouguer, mas deixa larga margem ao arbítrio 
de quem o escolhe. 

Recentemente, Oldham e Sutherland [34], para 
fazerem o estudo de um extenso campo gravítico 
regional, no Sudoeste da província de Ontario, 
introduzindo a hipótese de que a anomalia local 
se comportava como uma variável aleatória do 
valor médio nulo com variância constante para 
todas as observações e utilizando o método dos 
mínimos quadrados, estabeleceram um critério 
para determinar o grau do polinómio interpola- 
dor do campo regional. 

Do ponto de vista teórico o processo é cati- 
vante, mas conduz a cálculos muito extensos, de 
modo que a sua utilização fica condicionada ao 
emprego de calculadores electrónicos. Há con- 
tudo que frizar que a extensão da área onde se 
aplicar este método deve ter certa influência na 
legitimidade das hipóteses introduzidas sobre o 
campo local. Quere dizer, pensamos que devem 
utilizar-se, nos cálculos, valores do campo original 
colhidos em uma zona tão extensa, que, de facto 
se possa considerar nela o campo local de valor, 
médio nulo e variância constante (1), 


CAPÍTULO IV 
À interpretação dos dados gravimétricos 


1— Definições e objectivos 


Com a interpretação pretende-se caracterizar 
as ocorrências do sub-solo responsáveis pelas 
anomalias gravíticas, quanto à localização, den- 
sidade, profundidade, forma e dimensões. 

Já referimos (Capítulo II) que existe uma inde- 
terminação intrínseca nos problemas da interpre- 


(1 Em sentido análogo, embora usando outros algo- 
ritmos, tinha já estabelecido Grant [36], em 1954, um 
método estatístico para determinar o campo regional. 
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tação e, por isso mesmo, no caso geral a carac- 
terização de todos aqueles elementos não é 
possível. Mas, muitas vezes, tenta realizar-se 
a interpretação porque, além do conhecimento 
das anomalias gravíticas, se dispõe de certos 
dados fornecidos pela Geologia e outros métodos 
geofísicos, pelas sondagens ou por trabalhos 
mineiros, dados que limitam as configurações 
possíveis a alguns casos pouco numerosos. Este 
conhecimento estranho aos dados gravíticos, pode 
consistir em: 


1. Informação sobre a profundidade provável 
da configuração potenciante; 


2. Conhecimento da sua densidade; 


3. Informação sobre a forma provável da 
mesma. 


Não sendo conhecidas, nem a profundidade, nem 
a densidade, nem a forma, a interpretação limi- 
ta-se a localizar em planta anomalias densitárias 
e, quando muito a dar-lhe o sinal do contraste. 
E não se julgue que, mesmo em casos como este, 
a informação geofísica seja necessariamente de 
somenos importância. Antes pelo contrário, na 
grande maioria dos trabalhos de prospecção do 
petróleo e dos minérios a simples localização 
de uma anomalia densitária com o respectivo 
sinal do contraste, numa região, onde a Geologia 
nada podia afirmar, é de importância enorme 
e pode ter valor económico considerável, por- 
quanto se traduz na definição de uma área res- 
trita, onde haverá que sondar com elevada pro- 
babilidade de encontrar algo (1). Pode mesmo 
afirmar-se que a finalidade da prospecção geofi- 
sica é localizar, em planta, as áreas de anomalias 
densitárias [28]. 

Contudo, se condições geológicas ou outras 
justificarem um pormenorizada e rigorosa pros- 
pecção geofísica por gravimetria e se se dispuzer 
de dados estranhos, nomeadamente a profundide 
e o contraste densitário, poderá haver interesse 
em realizar a interpretação, tentando traduzir as 
anomalias por prováveis configurações do sub- 
-solo. 

Para tanto, existem diversas técnicas de inter- 
pretação. 


(1) Na prospecção do petróleo, os furos dirigidos pela 
geofísica, em 25 "/y dos casos são produtivos /28). 


2 — Técnicas de interpretação 


Uma vez que os elementos de campo sejam 
simultâneamente pormenorizados e precisos, pode 
tentar-se a interpretação seguindo dois métodos. 

Num deles admite-se que a configuração poten- 
ciante tem determinada forma, certo contraste 
densitário e se encontra a dada profundidade e 
calcula-se a anomalia correspondente, que se 
compara, em seguida, com aquela que foi obser- 
vada. Se a concordância é boa, poder-se-á afir- 
mar que a configuração hipotética é, tanto mais 
provável, quanto maior for a segurança das 
hipóteses que levaram à sua construção. Se a 
concordância não é boa, há que alterar alguns 
ou todos os elementos de configuração poten- 
ciante hipotética, calcular novamente o seu campo 
e tornar a fazer a comparação com a anomalia 
observada. E assim sucessivamente. 

Este método denomina-se de interpretação indi- 
recta. 

No outro método, chamado de interpretação 
directa, admite-se que a configuração potenciante 
tem determinado densidade ou pertence a certo 
tipo (anticlinal, massa, filão, falha, camada hori- 
zontal ou inclinada, etc.) e que jaz a certa pro- 
fundidade e, por intermédio de técnicas analíti- 
cas especiais, determina-se a sua forma e dimen- 
sões. 

Em qualquer dos dois casos a interpretação 
indicará a configuração potenciante mais prová- 
vel, em face dos dados gravíticos ou doutra natu- 
reza de que se dispõe. 


2.1 — Interpretação indirecta 


Pondo de parte o método de interpretação 
indirecta que consiste em supor a configuração po- 
tenciante com formas geométricas simples (cilin- 
dro horizontal infinito, cilindro vertical, esfera, 
estrato horizontal semi-infinito, etc.) os campos 
gravíticos das quais se podem deduzir por via 
analítica, mas cuja semelhança com as ocorrên- 
cias reais é, quase sempre, muito precária, pode 
dizer-se que o método mais empregado é o que 
utiliza ábacos, para calcular o campo devido à 
configuração hipotética. 

Podem construir-se ábacos próprios para as 
chamadas configurações planas ou bidimensio- 
nais, isto é, estruturas em que uma dimensão 
excede consideravelmente as outras, como as 
dobras de charneira muito longa, as falhas, os 


filões de grande extensão, as massas muito alon- 
gadas, etc., e que, portanto, provocam anomalias 
de características constantes segundo uma di- 
recção. 

Nestes casos, desenha-se uma secção hipoté- 
tica normal à maior dimensão da configuração e 
sobrepõe-se-lhe um ábaco traçado em papel trans- 
parente, fazendo coincidir a origem do mesmo com 
os pontos onde se pretende conhecer a anomalia. 

Estes ábacos são construídos, dividindo uma 
área, normalmente semi-circular, em sectores de 
dimensões tais, que a componente vertical do 
campo que cada um deles produz no centro do 
semi-circulo (origem do ábaco) é a mesma. Por 
consequência, fazendo a contagem do número de 
sectores que caem dentro da secção da configu- 
ração hipotética, obtém-se a gravidade no ponto 
coincidente com a origem do ábaco, a menos um 
factor, . que depende da escala em que se dese- 
nhou aquela secção, 

Quando a configuração potenciante não pode 
ser considerada bidimensional, há que realizar 
certas correcções ou utilizar ábacos para tal pre- 
parados [5], ou então dividi-la por planos hori- 
zontais, em «estratos» de altura constante, nos 
quais se supõe distribuída uma densidade super- 
ficial igual ao produto do contraste densitário 
admitido pela altura dos estratos e utilizar ábacos 
apropriados, que fornecem a gravidade provocada 
por qualquer dos estratos. Há aqui, contudo, que 
variar a escala da representação de cada estrato, 
porquanto um parâmetro que intervém no dese- 
nho do ábaco é a profundidade, facto que pode 
tornar o método trabalhoso. 

Com as configurações bidimensionais, usam-se 
também, por vezes, integradores mecânicos para 
determinação do campo. 

Além destes processos geométricos e mecãá- 
nicos, que referimos por alto, parece-nos de par- 
ticular interesse um outro baseado, em uma ana- 
logia óptica, cujo estudo tivemos ocasião de 
continuar durante a elaboração deste trabalho e 
ao qual fazemos referência mais detalhada no 
Capítulo VI. 

Com efeito, pode-se por: 


 Q 1 
dc dendçs Der (1), dd 
ou 
cos (z, 1) 
g (P) = k | d 
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Dividindo portanto, a configuração potenciante 
hipotética em estratos horizontais de altura dt 
fisicamente pequena, aos quais atribuímos a den- 
sidade de volume y, cada estrato pode ser asse- 
melhado a uma distribuição de simples camada 
com 7==y Àt e teremos para o estrato i: 


g, (P) =k | y At. sos (e, 1) dS (1-IV) 
i 


Mas o integral de (1-IV) representa, a menos 
um factor, a intensidade de iluminação em P 
provacada pela superfície S;, onde existe distri- 
buíida determinada intensidade luminosa (ver 
Apêndice). 

Esta analogia foi, parece-nos, pela primeira 
vez apontada por Tsuboi e Fuchida em 1938, 
que, tendo colocado uma chapa fotográfica a 
certa distância e paralelamente a uma superfície 
luminosa em forma de L obtiveram na chapa 
manchas que se assemelhavam notávelmente aos 
contornos da anomalia gravítica calculada para 
estrato com aquela mesma forma. 

Quere dizer, se, em vez de utilizar uma chapa 
fotográfica se medir a intensidade de iluminação 
provocada num determinado plano por uma su- 
perfície plana luminosa paralela àquele e com a 
forma de qualquer dos estratos em que se divi- 
diu a configuração hipotética, obter-se-á, a me- 
nos uma constante, a distribuição da anomalia 
gravitica produzida pelo estrado considerado. 
Repetindo a operação para os estratos em que 
se dividiu a configuração hipotética, e somando 
em cada ponto os valores obtidos, construiremos 
de modo muito simples e também rigoroso a 
anomalia gravítica daquela. 

Este processo não envolve outra aproximação 
que não seja a resultante da divisão da massa 
potenciante em estratos de altura fisicamente 
pequena, é susceptível de grande rigor, desde 
que se disponha de dispositivo sensível para 
medir intensidades de iluminação, como é o caso 
de um bom fotómetro, e pode utilizar-se com 
qualquer tipo de configuração potenciante. 

É além disso fácilmente adaptável a um braço 
mecânico que se desloque sobre uma carta tra- 
çada em qualquer escala, onde estejam marcados 
os pontos nos quais se deseja conhecer o valor 
da anomalia, de modo que a leitura das intensi- 
dades luminosas, para os diferentes extractos se 
faça sempre nos mesmos pontos. 
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Em resumo, trata-se de um processo labora- 
torial para o estudo do campo gravítico, utili- 
zando um modelo óptico barato, onde a relação 
de semelhança é muito simples, e em que o tra- 
balho dispendido para a estruturação do campo 
é bastante reduzido. Contamos, por isso, poder 
desenvolvê-lo e aperfeiçoá-lo. 


2.2 — Interpretação directa 


Como referimos, com a interpretação directa 
pretende-se, partindo do conhecimento da ano- 
malia, estabelecer a configuração potenciante 
anómala. Vimos também que só é de tentar 
esta interpretação, quando se dispõe de dados 
gravimétricos abundantes e precisos e ainda, se 
se puderem estabelecer certas hipóteses sobre a 
profundidade, o contraste densitário e o tipo de 
configuração que se procura. 

A interpretação directa realiza-se por meio do 
«prolongamento descendente» (downward conti- 
nuation). 


2.2.1 — O prolongamento descendente 


O prolongamento descendente consiste em 
calcular o campo tão próximo quanto possível 
da massa potenciante, uma vez conhecida a ano- 
malia na superfície da Terra. 

A razão de ser deste modo de proceder ba- 
seia-se nas seguintes considerações: 

Vimos, no Capítulo II, que sobre a camada 
equivalente, os valores da densidade superficial 
estão relacionados com a gravidade por: 


as IR) 
2 mk 


onde g (Q) é o campo, na camada equivalente, 
e k a constante de gravitação. 

Mostramos também que aquela camada podia 
ser substituída por uma estrutura equivalente 
com um único contraste desintário, desde que, 
se fixasse este em certo valor 7, e que o relevo 
dessa estrutura, referido ao plano da camada 
equivalente seria então: 


À gs MO 
2 TT kyj 


Isto é, a menos um factor contante, o campo 
na vizinhança da estrutura equivalente, é igual 


